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Abstract 

We study the low-temperature properties of S = 1 and 1/2 alternating spin 
chains with antiferromagnetic nearest-neighbor exchange couplings using an- 
alytical techniques as well as a quantum Monte Carlo method. The spin-wave 
approach predicts two different low- lying excitations, which are gapped and 
gapless, respectively. The structure of low-lying levels is also discussed by 
perturbation theory in the strength of the Ising anisotropy. These analytical 
findings are compared with the results of quantum Monte Carlo calculations 
and it turns out that spin-wave theory well describes the present system. We 
conclude that the quantum ferrimagnetic chain exhibits both ferromagnetic 
and antiferromagnetic aspects. 
PACS numbers: 75.10Jm, 75.30Ds, 75.40Mg 
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I. INTRODUCTION 



The low temperature properties of low-dimensional quantum antiferromagnets have been 
of great interest for many years, in particular since Haldane [ 1] made the prediction that 
integer-spin and half-odd-integer-spin Heisenberg antiferromagnetic chains should behave 
very differently. Now, with the help of analytical methods and various numerical approaches, 
it is well accepted that the integer-spin chain is massive, whereas the half-odd-integer- 
spin chain is massless. This stimulated also several attempts to investigate the quantum 
behaviour of chains consisting of two types of spins. An integrable model of this type 
was constructed by de Vega and Woynarovich [ §], which allows us to guess the essential 
consequences of chains with spins of different length S. Recently, several authors [ ||],|4| 
discussed in detail such a chain with spins S = 1 and 1/2 as the simplest case. However, these 
integrable models usually include complicated interactions and very little is known about 
the pure Heisenberg model which is supposed to describe real ferrimagnetic compounds [ [5] . 

In the following we therefore study the S = 1 and 1/2 Heisenberg antiferromagnetic 
chain defined by the hamiltonian (|l|) 

L 

n = j^[(s l -s l ) x + (s l -s l+1 ) x \ , (i) 
i=i 

where (Si ■ Sj) x = A( Sf s] + SJs) ) + Sf s] with Si and S; being the S = 1 and S = 1/2 
spin operators, respectively, and L is the number of unit cells. We adopt periodic boundary 
conditions and investigate ground-state properties as well as the low-lying excited states 
spin-wave theory as well as series expansions with regard to the Ising anisotropy A and a 
quantum Monte Carlo (QMC) method. 

Basic predictions on the low lying level structure can be made by applying the Lieb- 
Mattis- Theorem [ ^| and the Goldstone-Theorem [ [F]. First we consider the ordering of the 
energy levels. The A sublattice is defined by the S = 1 sites (maximum spin is Sa = L), the 
B sublattice is defined by the spin-1/2 sites (maximum spin is Sa = L). The Hamiltonian 
([[]) consists of intersublattice interactions only, which allows us to apply [ || to our system. 
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The energy levels order in the following manner: 

E(S + 1) > E{S) for all S > S E{S) > E{S) for S < S 

Here S = \Sa — Sb\ denotes the total spin of the ground state which takes the value 
S = L — L/2 = L/2 in our case. Therefore the ground-state degeneracy is L + 1. Although 
the ground states have finite spin S the following points are in contrast to the usually 
ferro magnet: 

• The ground-state degeneracy of the ferromagnet is larger by a macroscopic amount. 

• Due to quantum fluctuations the ground state deviates from the ferrimagnetic Neel 
state, | + 1,-1/2, +1,-1/2,...). 

Because each of the ground states of the chain breaks the rotational symmetry of the Hamil- 
tonian, we can apply the Goldstone-Theorem [ 0] to predict a gapless excitation. For mag- 
netizations lower than L/2 the ferrimagnet is comparable to a ferromagnet and therefore 
this gapless excitation should belong to a branch of excitations of ferromagnetic character. 

The paper is organized as follows: In the second section we present the spin-wave ap- 
proach to calculate the dispersion relations as well as the ground-state energy and the 
ground-state correlation functions. In the following section we study the low-temperature 
properties employing a QMC method. The numerical data will be compared to the results 
of the first section and to perturbation calculations. Conclusions will be given in the final 
section. 

II. SPIN WAVE THEORY 

In this section the spin wave theory (see e.g. [ ||) is applied to our system in the large 
S limit. For this purpose the odd sites are assigned spins gS and the even sites spins S. 
To discuss the case of our Hamiltonian (|I]) we use g = 2 large S limit. We start from the 
fully ordered state with Sf ot = M = L/2 and use the following spin operators in the two 
sublattices in the lowest order in 1/5 : 
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sublattice A : S* = gS — a} n a n , = \J2gSa n , (2) 

sublattice B : s z n = -S + b% s+ = V2Sb{ , (3) 

Here a' n , b^ n are Bose operators. We expect the spin wave theory to give qualitatively 
correct results for the following reason: In the classical antiferromagnet the two Neel states 
become disconnected by the transformation (|2],|3|) since the action of the transverse Hamil- 
tonian on one of the two Neel states does not lead to the other Neel state. However, it 
is known that domain wall excitations, namely excitations connecting the two Neel states, 
are very important for the spin- 1/2 antiferromagnet. Spin wave theory is wrong for the 
spin-1/2 antiferromagnet since domain wall excitations cannot be taken into account. In 
the spin 1-1/2 system the two Neel states are disconnected in principle because they belong 
to subspaces with different magnetization. Thus a classical ground state including quantum 
fluctuations should be qualitatively correct. 

In the limit S — > oo the interaction terms in the Bose Hamiltonian become negligible and 
we end up with the following expression for the Hamiltonian, bilinear in the Bose operators, 

Hsw = ~2gS 2 JL + 2SJ]T [gb\b k + a\a k + ^cos{k){a\b\ + a k b k )\ , (4) 

k 

where a) n = L -1 / 2 J2 k e lkn a\. and tf n = L -1 / 2 J2 k e~ lkn b\. The Hamiltonian is straightforwardly 
diagonalized into 

Hswt = E + 2SJ ajafc + ^tAPi) > ( 5 ) 

k 

= \ (j(g + iy-4gcos*{k) ±(g- 1)) , (6) 



E = JSJ2 



[l + gy-4gcos*(k)-(l + g) 



gJS 2 N, (7) 



where the eigenvectors are determined by 



a\ = cosh(%)4 + sinh(?7 fc )6 fe , p\ = smh(r] k )a k + cosh(r] k )bl , (8) 
tanh(2%) = cos(A;) . (9) 

The results of the spin wave theory contain both ferromagnetic and antiferromagnetic as- 
pects: In the ferromagnetic branch (u^) we obtain a gapless spin wave with M — L/2 — 1. 



For small values of the wave vector the dispersion is uo^ — g/(g — 1) k 2 + . . . The quadratic 
behaviour with wave vector k indicates the ferromagnetic character of this mode. 
The antiferromagnetic spin wave with M = L/2 + lis gapped. The magnitude of the gap is 
exactly J when we put g — 2, S — \. This result will be compared to the result of the QMC 
calculation to be presented in the next section. If we chose the spins such that the spin 
lengths in the different sublattices are equal, i.e. g — > 1 the gap of the antiferromagnetic 
branch vanishes. Another interesting feature is the ground-state sublattice magnetization 
as defined by 

L 

(M A ) g . s . = L - £ 4«n = L(l - r) (10) 

n=l 

t = - dksmh 2 (ri k ) « 0.305 . (11) 

7T JO 

Notice that the one dimensional ferrimagnet has a finite spin reduction r. This is due to 
the factor 2y / 2/3 < 1, obtained for g = 2 in relation (Q. In the usual antiferromagnet the 
corresponding factor is 1 and leads to a diverging spin reduction (in 1 D). 
The ground-state energy per unit cell in the spin wave approach is obtained as 

E /J = ~L + — r dkJ9 - 8 cos 2 (A;) w -1.4365L (S = \,g = 2). (12) 
2 2,71 Jo 2 

Further we studied the ground-state correlation functions between two sites with S = 1 

for the ground state with magnetization L/2 and obtained the following results: 

(S*S* +r ) = (1 - r) 2 + f(r) (13) 
(S+S~ +r ) oc (14) 
f( r ) = J2 E ei2 " r smh 2 ( % ) cosh 2 ^,) (15) 

k,q 

The asymptotic behaviour for the function f(r) is calculated by taking the continuum limit 
for the sum in eq.(^5|). We observe asymptotical exponential decay with a correlation length 
£, = 21n(2). So for the intrasublattice distance r > 2 correlations have decayed and the 
square of the spin reduction (1 — r) 2 remains. From comparison with the QMC results it 
will be shown later this is a qualitatively correct picture. 
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As conclusion of this section we discuss the case of strong alternation, a limit which 
naturally reproduces the scenario as given above. If we introduce the interaction strength 
J between sites (2n) and (2n + 1) and 5 J between sites (2n — 1) and (2n) the Lieb-Mattis- 
Theorem and the Goldstone-Theorem still hold. The M = L/2 ground state in the dimerized 
limit 5 = can be written down as 

\0) = \D + ) 1 \D + ) 2 ... \D + ) L 
(S n + s n f\D + ) n = \ Q + l) \D + ) n {S z n + s z n )\D + ) n = l -\D + ) n 

\D+) = ^(10, +1/2) -y/2\ + 1,-1/2) 

For values <5 <C 1 the system behaves like a S = \ ferromagnet with the dublets (D + ,D_) 
being the effective spin S e ff — \. A first order perturbation calculation in the alternation 
parameter 5 leads to a ground-state energy of E = — L(l+l/9)+0(5 2 ). Now a ferromagnetic 
spin wave can be constructed within the dublet subspace 

l^^fe^K.-.m-PA (16) 

The dispersion up to the first order in 5 is uo(q) = 4/9A(l — cos(g)) + 0(5 2 ). It is gapless 
and proportional to q 2 for small wavevectors as one would expect. 

In addition to the ferromagnetic excitations we can construct antiferromagnetic spin waves, 
i.e. spin waves with S^ ot = L/2 + 1. Again we discuss the case of strong alternation and end 
up with 

\q) = ^=i:e^\D + ) 1 ...\Q 3 ) n ...\D + ) L (17) 

n =l 

Here \Q 3 ) denotes the quartet with magnetization 3/2 . The first order dispersion is 
u(q) = 3/2 + A(7/18 - 2/3cos(g)) + 0(5 2 ) and is gapped. 



III. NUMERICAL RESULTS 
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A. Brief Account on the Numerical Procedure 



In the following we employ a quantum Monte Carlo method based on the Suzuki- Trotter 
decomposition [ Q of checker-board type [ 10 1. Raw data are taken for a set of Trotter 
numbers n and are extrapolated to the n — > oo limit with the parabolic fitting formula. We 
carry out all the calculations in certain subspaces with a fixed value of the total magneti- 
zation. Since we treat the chains with periodic boundary conditions, not only the Monte 
Carlo flips of local type but also the global flips along the chain direction are taken into 
the numerical procedure. On the other hand, global flips along the Trotter direction, which 
mean fluctuations of the total magnetization, are not included in order to well describe the 
ground-state properties. The quantum Monte Carlo algorithm to update the spin config- 
uration is detailed elsewhere [ [IT[| . We have confirmed that almost the same results are 
obtained at two temperatures, k-^T/J = 0.04, 0.02, and thus we regard these temperatures 
as low enough to successfully extract the lowest-energy-state properties. Here we show the 
data taken at k-^T/J = 0.02. The data precision is almost four digits for the energy, and 
two digits for the spin correlations. 



B. Low-Energy Structure 

In Fig.[I] we plot the lowest energies per unit cell in the subspaces with M = J2i(S- +s^) = 
L/2, L/2 — 1 as a function of L. The coincidence between both sets of data is nothing but 
the numerical demonstration of the above mentioned Lieb-Mattis theorem, that is, the L/2- 
multiplet structure of the ground state. The rapid convergence into the long-chain limit 
suggests a rather small correlation length in this system, which will be actually observed 
in the following. Within the present numerical precision, the ground-state energy in the 
thermodynamic limit is estimated as Eq = —1.455 ± 0.001. As we have observed in Fig.[l|, 
even the ground-state energy of the L = 16 chain is already close to the thermodynamic-limit 
value. Thus we have plotted the quantum Monte Carlo data for L = 16. 
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Although the quantum Monte Carlo data as presented above are already conclusive in 
themselves, we present an additional argument based on perturbation theory in the Ising 
anisotropy A. We compare in Fig.0 the quantum Monte-Carlo estimated ground-state ener- 
gies and the corresponding perturbation-theory result 

\2 \4 

- Eq/L = 1 + — - — - 0.05136A 6 + 0.02809A 8 + C(A 10 ) , (18) 
2 48 

as a function of A. We find a fairly good agreement between the eighth-order perturbation 
result and the QMC calculation. What should be emphasized is that fourth order calculation 
shows good agreement with the correct result. This fact gives us an idea on the spin 
configuration in the M = L/2 ground state. All the fluctuations introduced within the 
fourth order are essentially classified into the following three types: 

(a) 2-sites fluctuation : 1, -1/2, 1, -1/2, 0, 1/2 , 1, -1/2, 1, -1/2, • • ■ 

(b) 3-sites fluctuation : 1, -1/2, 1, 1/2, -1, 1/2 , 1, -1/2, 1, -1/2, • • ■ 

(c) 4-sites fluctuation : 1,-1/2, 0, 1/2, 0, 1/2 , 1, -1/2, 1, -1/2, ■ • - 

We note that the formation energy of the defect (c) is not twice as much as that of (a) and 
therefore they should be distinguished. Based on these fluctuations we are led to discuss 
microscopic quantum fluctuations at the isotropic point (A = 1). We present in Fig.|| a QMC 
snapshot from which we can extract the image of quantum fluctuations to a certain extent. 
Here we show the snapshot at k^T / J = 0.02 for A = 1.0, L = 32, and n = 20, where the 
horizontal and vertical lines denote the chain and the Trotter directions corresponding to 
space and time, respectively. We find out everywhere local defects breaking the Neel Order, 
whereas they are all identified with the above-mentioned fluctuations. Thus we expect the 
fourth order calculation to well describe the the M = L/2 ground state. 

Next we discuss the lowest excited states for the ferromagnetic and the antiferromagnetic 
branches. In the subspace with M = L/2 — lwe construct a magnon for A = by flipping 
a spin-1: S z = 1 — > S z = 0. Expanding this state and taking the limit k — > we obtain 
the lowest mode in this subspace. Such an expansion up to the fourth order is compared 
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with Eq{L/2 — 1) — Eg{L/2) from the QMC calculation in Fig.[|. The perturbation result 
gives reasonable results up to A ~ 0.8. A comparison between them suggests that the k = 
magnon condenses into the ground state with M = L/2 — 1 at the isotropic point. 
The antiferromagnetic magnon is constructed by flipping a spin- 1/2 up. Again the pertur- 
bation result and the QMC data for Eq{L/2 + 1) — Eq(L/2) are illustrated in Fig.[|. Here 
the validity for the fourth order expansion only reaches up to A ~ 0.4. From the QMC 
calculation we obtain a gap of A/ J = 1.767 ± 0.003. Here somewhat larger uncertainty, 
rather than one for the ground-state energy, mainly comes from the Monte Carlo estimate 
in the subspace ofM = L/2 + l. This value is much bigger than the prediction of the spin 
wave theory A = J. 

We conclude this section by pointing out that the quantum behaviour of the 1-1/2 system 
results in an enhancement of the gap. From Fig]4] we observe that the pure Ising energy of 
2 J is just lowered by a small amount when moving to the isotropic point. 

C. Spin Correlations 

We show in Fig.|5| the spin correlation functions between the spins of the same type for 
the lowest-energy states in the subspaces with M = L/2 — 1, L/2, L/2 + 1. We note that the 
self correlation of spin-1 deviates from 2/3 because of the multiplet structure of the ground 
state. In comparison with the spin correlations for M = L/2, those for M = L/2 — 1 are 
significantly reduced rather than those for M = L/2 + lin the case of spin-1, and vice versa 
in the case of spin- 1/2. This is well understood considering that the ferromagnetic and the 
antiferromagnetic magnons, which exist in the subspaces with M = L/2 — 1 and M = L/2+1, 
originate from the spin flips in the S — 1 and the S = | sublattices, respectively. Those 
excitations are expected to reduce the ferromagnetic correlations between spins with S = 1 
and spins with S = \. We note that in the thermodynamic limit, both spin correlations for 
M = L/2 — I and M = L/2 + 1 should coincide with those for M = L/2. Nevertheless, Fig.| 
is still useful because it suggest to a certain extent the thermodynamic-limit spin correlations 
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in the subspaces with magnons of finite density. 

In the above sense, now let us concentrate ourselves on the subspace with L/2. We 
have already confirmed that the spin wave theory gives the asymptotic exponential decay 
of the spin correlations, where the correlation length is estimated to be less than unity. We 
here observe so rapid decay of the correlations that an estimate of the correlation length 
is beyond the present numerical precision. However, careful observation of Fig.|| shows us 
that the correlations between spins S — ^ are a little bit less rapid than those between spins 
S — 1. Spin wave theory cannot reproduce this feature because the decay is determined 
by the function f(r) in eq . fll5|) for both sublattices. We further point out that the Monte 
Carlo calculation gives the spin reduction r as r ~ 0.21, which is somewhat smaller than 
the spin- wave-theory result. Therefore the spin- wave-theory cannot quantitatively describe 
the quantum fluctuations. 

IV. CONCLUSION 

We have calculated ground-state properties and low-lying excited states for an alternating 
ferrimagnetic spin chain with spins S = 1 and ~. The ground state is a spin S = L/2 multi- 
plet. The model consist of a ferromagnetic and an antiferromagnetic branch corresponding to 
magnitizations M < L/2 and M > L/2, respectivly. The ferromagnetic branch has gapless 
excitations with dispersions uj oc k 2 , k <C ir according to spin wave theory. The antiferro- 
magnetic branch with M = L/2 + 1 shows a gapped spin wave with A/J= 1.767 ± 0.003. 
Both branches have longitudinal correlation functions consisting of a constant (square of 
spin reduction) plus strong exponential decay and therefore both branches show long range 
order. The manifestation of quantum behaviour shows up in the following points 

• the gap for the M = L/2 + 1 excitation is enhanced compared to the spin wave theory. 

• the L/2 ground state deviates from the Neel state owing to quantum fluctuations. 

• there is an indication for the correlation length to be larger in the S — 1/2 than S = 1 
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sublattice. 

In order to study the mechanism for the gap formation in detail we have to investigate the 
perturbation theory to higher order. This is under investigation as well as the construction 
of matrix product states as variational ground state for M = L/2 [ |T^]. The construction 
of these matrices can be based on the fact that unit cells with magnetization —3/2 show up 
rather rarely in Fig.^. Therefore the system exhibits a tendency to weak ferromagnetism 
similar to the one discussed by Niggemann and Zittartz [ |l3j in S = 3/2 chains with matrix 
product ground states. 
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FIGURES 

FIG. 1. Size dependences of the ground-state energy in the subspaces with M = L/2 and 
M = L/2 — 1 

FIG. 2. Dependences of the ground-state energy on the Ising anisotropy A obtained by a QMC 
method (O), and 4-th-order (dotted line) and 8-th-order (solid line) perturbation theory. 

FIG. 3. A snapshot of the transformed two-dimensional Ising system, where the horizontal and 
the vertical lines denote the chain and the Trotter directions corresponding to space and time, 
respectivly, and '+', '0', and '=' denote the spin projections +1, +1/2, 0, —1/2, —1 . 

FIG. 4. Dependences of the excitation gaps Eq(L/2 — 1) — Eq(L/2) (O) an d 
Eg(L/2 + 1) — Eq(L/2) (O) on the Ising anisotropy obtained by a QMC method, where the 
results within the 4-th-order perturbation treatment are also shown (solid lines). 

FIG. 5. The groud-state spin correlations between spins S = 1 (a) and between spins S = 1/2 
(b) in the subspaces of M = L/2, L/2 - 1 and L/2 + 1 for L = 16. 
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